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We examine the semiclassical content of SU{3) Yang Mills theory on the lattice at hnite tem¬ 
perature. Employing the cooling method, a set of classical helds is generated from a Monte Carlo 
ensemble. Various operators are used to inspect this set with respect to topological properties. We 
hnd pseudoparticle fields, so-called caloron solutions, possessing the remarkable features of (super¬ 
positions of) Kraan-van Baal solutions, i.e. extensions of Harrington-Shepard calorons to generic 
values of the holonomy. 


I. INTRODUCTION 

Guided by the idea that the semiclassical approach to 
QCD 0,0) 13 provides a reasonable model for the QCD 
vacuum, one first has to specify the basic set of classical 
solutions of the underlying Yang-Mills theory. Mathe¬ 
matical rich and interesting structures naturally arise in 
the construction and topological classification of these 
fields. For SU{N) gauge theory in the Euclidean space 
R'* solutions with finite action are the well known BPST 
instantons 0 and their generalizations to arbitrary topo¬ 
logical charge [3 . The vigorous work pursuing the semi¬ 
classical approach culminated in the so-called instanton- 
liquid model 0 0 13 > which successfully realizes funda¬ 
mental non-perturbative properties of QCD such as the 
U{1)a anomaly and chiral symmetry breaking. It has led 
to numerous phenomenological applications (for newer 
reviews see 11113 ). However, it failed to reproduce one 
essential feature of QCD - confinement.^ 

In case of the SU(N) gauge theory with finite temper¬ 
ature defined on x quark (de)confinement is closely 
related to the global center symmetry. The Polyakov loop 
at temperature T 

/ ,h=\/T 

P(x)=Pexp / AQ{xo,x)dxo 

an order parameter of this symmetry, signals a phase 
transition from a deconfining to a confining phase. De¬ 
confinement is a feature of the high temperature phase, 
where the symmetry is spontaneously broken and ~ in 
the limit of large temperatures - the Polyakov loop takes 
values in the center of the gauge group Z{N). Below 
the critical temperature the symmetry becomes restored, 
the trace of the Polyakov loop fluctuates around zero and 
confinement is observed. One can hope that not only for 
temperatures above the transition temperature the semi¬ 
classical approach is justified. Then it should be able to 


^ This incited the work of Negele et ah, who were able to repro duc e 
confinement in a gas of regular gauge instantons or merons 


reflect the transition as an important feature of the finite 
temperature gauge theory. 

A semiclassical model for non-zero temperature 
QCD has first been formulated in m starting from 
Harrington-Shepard (HS) caloron solutions which 
represent infinite chains of time-periodic instanton solu¬ 
tions in R'^. Restricted to the strip x R^ HS calorons 
have an integer-valued topological charge Qt G Z like 
BPST instantons. For SU{N) HS calorons are realized 
as SU (2) embeddings. More restrictive from our point 
of view is the fact that their vector potentials behave 
asymptotically such that trivial holonomy emerges, i.e. 
their asymptotic Polyakov loop values belong to the cen¬ 
ter of the group 

iPoo = Jim r{x) G Z{N). (2) 

|x|—»^oo 

If one would have to describe only the deconfinement 
phase, this asymptotic behavior which is directly related 
to the broken center symmetry values of the order pa¬ 
rameter is certainly adequate. Is this the case also for 
the confinement phase ? The authors of have ar¬ 
gued that the path integral evaluated in the infinite vol¬ 
ume limit cannot be dominated by classical fields hav¬ 
ing a non-trivial holonomy. But this conclusion does not 
hold for a single topological excitation created in a fi¬ 
nite box or in the case of a finite density of such ob¬ 
jects. Therefore, the question arises, whether there are 
also classical fields asymptotically matching to a van¬ 
ishing order parameter as required in the confinement 
phase. Indeed, selfdual time-periodic caloron solutions 
of the SU{N) Yang-Mills equations with arbitrary (non¬ 
trivial) holonomy have been explicitly constructed by 
Kraan and van Baal M , and Lee and Lu M several 
years ago. The new SU(N) calorons - we call them KvB 
calorons - with integer-valued topological charge Qt 

consist of N IQJ _ ^constituents, which appear as static 

BPS monopoles [l^ Il7| if sufficiently separated. Their 
masses T izm are then determined by the differences 
of holonomy eigenvalues Vm = Pm+i — Pm , 

1^00=9 exp(27ri diag(/xi,/r 2 ,..., pn)) 5^ (3) 

ordered such that pi < p 2 < • ■ • < Pn+i = ^ + Pi 
and + ^2 + ■ ■ ■ + Pn = 0- The separated and massive 
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caloron constituents can be localized by one or more of 
the following criteria [T^: 

(a) at the centers of their three-dimensional spherical 
lumps of action or topological charge, 

(b) at the positions, where the zero-mode density of 

the Dirac operator in the background of the caloron 
field Ip = 7^(9^ -I- is maximal, depending 

on the temporal boundary condition for the fermion 
field, 

(c) at the points x, where at least two of the eigenvalues 
oiV{x) coincide, 

(d) at the positions of static Abelian (anti)monopoles 
representing defects of an appropriately chosen 
Abelian gauge. 

However, if the constituents are close to each other or 
if they become even massless, these criteria do not hold 
unambiguously. The first two criteria cannot distinguish 
between the constituents at all. The constituents form 
joint non-static caloron lumps of topological charge. The 
last two criteria are still working, but the extracted po¬ 
sitions will not completely agree. 

In any case, the non-trivial holonomy suggests a proper 
link between the classical background field and the dy¬ 
namics of the order parameter. Therefore, a semiclassi- 
cal (variational) approach to Yang-Mills theory should be 
based on the new, KvB calorons. That such an approach 
is promising has been recently argued by Diakonov M- 
He and his collaborators have shown that the quantum 
amplitude in the background of KvB calorons leads to 
interactions of the monopole constituents at the quan¬ 
tum level Ellis These interactions become attractive 
at higher temperature (i.e. deconfinement) and repul¬ 
sive at lower temperature (confinement). Thus, in the 
deconfinement phase HS calorons might be favored (with 
all but one constituent massless, i.e. invisible), whereas 
in the confinement case the repulsive interactions sup¬ 
port a view of a gas to be described fully in terms of 
monopoles or KvB caloron constituents. The calculation 
of the fluctuation determinant still awaits extension from 
the SU{2) case El H to arbitrary N. The zero-mode 
contribution to the semiclassical integration measure can 
be obtained for arbitrary N from the determinant of the 
explicit moduli space metric as computed by Kraan a few 
years ago [l^. This, as well as a more detailed analysis of 
how in particular the parameters of the N monopole con¬ 
stituents are related to the scale parameter of the stan¬ 
dard SU{N) instanton, can be found in Ref. 

In recent papers El published together with 
Martemyanov, Shcheredin and Veselov we have explored 
the existence of KvB calorons in the (simpler) case of 
SU{2) lattice gauge theory. We first relied on the cooling 
method in order to search for (approximate) solutions of 
the Yang-Mills equations of motion. Later on we have 
investigated smoothed configurations obtained by four¬ 
dimensional smearing i.e. “keeping closer” to equilibrium 


fields El- In both cases, starting from Monte Carlo equi¬ 
librium fields generated within the confinement phase we 
have shown that non-trivial holonomy with tr'Poo ~ 0 
dominates the vacuum structure. The topological lumps 
we have found were typically composed of monopole con¬ 
stituents and possessed all properties of KvB solutions 
(or their constituents). 

In this paper we are extending our investigations to 
the SU (3) case of finite temperature lattice gauge fields 
in the confinement phase. Relying on the cooling pro¬ 
cedure our aim is again to explore the structure of 
(semi)classical fields at the level of a few instanton ac¬ 
tions S = n S'inst, n = with S-mst = 87r^/(/^. 

For this we had to refine our methods and to improve 
lattice observables. Our paper is organized as follows. 
In Section HD we review known analytical formulae for 
the SU{3) KvB caloron.^ Later on, in Section ITTTl we 
describe the observables and the methods which are 
used to compare the analytical results with the classi¬ 
cal lattice gauge fields. Results for typical examples of 
cooled lowest-action (classical) fields are presented in Sec¬ 
tion HU The results obtained for ensembles of classical 
fields extracted also at higher action level are discussed 
in Section El 

II. PROPERTIES OF KVB CALORONS 

Most analytical efforts are currently directed to im¬ 
proving the knowledge of higher charge SU{2) calorons 
[30j. Expressions are known for topological c harg e-one 
calorons in SU{2) EJ and for generic SU{N) We 
recall the expressions of the action .s{x) and the den¬ 
sity of the fermionic zero-mode ofthe5'C/(3) 

single-caloron field with jQtl = 1 ■ The zero-mode 
depends on the temporal boundary condition imposed, 
'I’zi.xo + b,x) = z{xo, x) (with b = 1/T). 

s(a;) = -]^d‘^d'^\nip{xo,x) (4) 

vl/tvl/,(a;) = -^aV.(A^) (5) 

+ ( 6 ) 

The auxiliary fields in Eqs. l|3Hni) are defined in the Ap¬ 
pendix. The caloron moduli space consists of three (con¬ 
stituent) monopole positions ijm in (m = 1,2,3 
and 2/4 = yi) and three phases in U{1) x [7(1) of which 
one plays the role of the overall time location. The 
temporal extension is scaled to b = 1 and the holon¬ 
omy is left arbitrary. For illustration the action density 
for a SU (3) KvB solution with well-dissociated monopole 
constituents is plotted in Fig.^ 


^ For further re ading about calorons consider the literature, e.g. 

Hi, HI HI m 
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FIG. 1: Logarithmic plot of the action density of a SU{3) 
caloron (from (s^)- 


zero-mode jumping. In the deconfined phase, they found 
a wide range in x where the zero-mode didn’t change 
localization, becoming delocalized only at one angle z 
corresponding to the phase of Loo- The localization of 
the zero-mode wave function was in agreement with lo¬ 
cal clusters of equal sign of the topological charge den¬ 
sity [S^. Globally, however, the latter was far more com¬ 
plex than that of a single caloron. 

Notwithstanding this state of affairs it is suggestive 
to carry over the lattice observables developed for the 
purpose of analyzing classical configurations to a pro¬ 
gram to rediscover superpositions of KvB calorons and 
anticalorons in Monte Carlo ensembles. Our aim is even¬ 
tually to prove the validity of the semiclassical approach 
near to the transition temperature and above. 

III. LATTICE TOOLS OF THE ANALYSIS 


Provided the constituents are well separated, the 
Polyakov loop values at their positions ijm, rn = 1,2,3 
are 


nm) 

= diag( 


g27r2/i3 

V{y2) 

= diag( 


g-TTZ/ii 

V{y3) 

= diag(-e-’"*^^ 



In Fig. 121 we illustrate the range of Polyakov loop values 
L{x) = (1/3) tTV{x) for two limiting cases: on the left 
hand side for a SU{3) KvB caloron with trivial holon- 
omy Loo = (1/3) tr^oo = +1 which corresponds to a 
SU{2) embedding of a Harrington-Shepard caloron ^ 
with real values of L(x) only, and on the right hand side 
for maximally non-trivial holonomy Loo = 0 leading to 
dissociated equal-mass constituents and L{x) populat¬ 
ing the inner triangle. 

The features of KvB calorons are the following. 


Working on the lattice has two obvious disadvantages 
- the explicit breaking of scale-invariance and the exis¬ 
tence of lattice artifacts. The classical action of a multi- 
caloron or instanton solution with topological charge Qt 
is S' = IQtl Sinst- It does not depend on the parame¬ 
ters of the solution: group orientations, positions and 
scale parameters. On the contrary, the standard Wilson 
action 

S'^ = ^ (1 - 1 Re tr t/,.^,), (8) 

depends on the scale size p of an instanton discretized 
on the lattice. It can be calculated as 113 


qW _ q 
■^inst ~ ‘Jinst 



Oia/pf 


(9) 


1. The holonomy can be chosen arbitrarily. 

2. The multi-caloron contains N jQtl constituent 
monopoles in accordance with various definitions. 

3. The fermionic zero-mode density 'I'll['^(a:) is local¬ 
ized at the constituent if z G {pm, Pm+i) and 
jumps (delocalizes) for z = pm- 

These remarkable features will be illustrated for calorons 
of various topological charge Qt, which are obtained on 
asymmetric lattices, i.e. Nr x with > Nr, by the 
cooling technique. This should be seen as an illustration 
of an extended set of caloron solutions. 

If realized on the lattice, the holonomy obviously can¬ 
not be chosen freely but dependent on the actual finite- 
temperature phase. Inspired by the properties listed 
above, Gattringer and Schaefer [Sa] have analyzed 
a subsample of Qt = ±1 Monte Garlo configurations, 
taken from the confinement and deconfinement phase, 
for the localization of the (single) zero-mode with chang¬ 
ing z . In the confinement phase, identified with max¬ 
imally non-trivial holonomy, they found the pattern of 


with d 2 = —1/5. Hence this lattice action favors small in- 
stantons. For a KvB caloron, the intrinsic scale is deter¬ 
mined by the constituent separations expressed by pm , 
and similarly small pm are preferred. For SU{2), in 
the limit of separated constituents, pm = d b/ir'^ with 
b = IjT and d being the distance between the two con¬ 
stituents. As a remedy one could choose an (improved) 
action with zero or even positive coefficient d 2 . For the 
purposes of this paper, however, we preferred to keep the 
numerically cheaper Wilson gauge action as the cooling 
action and simply to monitor the scale-size dependence. 

Due to the finiteness of the lattice one cannot repro¬ 
duce the X infinite toroidal or cylindric topology. 
Instead, depending on the aspect ratio Nr/Na one can 
interpolate between the topologies of the 4-torus T'^ and 
that of a cylinder. On the other hand, using the lattice 
as stage of such investigations has also advantages, for 
example the possibilities 

... to produce classical fields in different topological 
sectors without knowing explicitly the correspond¬ 
ing analytic expressions. 
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FIG. 2: Schematic scatterplots (left: on the real axis, right: inside the triangle) for the local Polyakov loop values L{x) in 
the complex plane, for a caloron with trivial holonomy L^a = +1 (left) and with maximally non-trivial holonomy Loo = 0 
(right). Black dots indicate the holonomy values, whereas white circles on the curved boundary give the Polyakov loop values 
at the monopole constituents according to Eq. 0 where two eigenvalues of V become degenerate in accordance with the 
monopole criterion (c) in Section 0 


... to include fermions into the study of solutions of 
the pure gauge theory, 

... to test the predictions of corresponding semiclassi- 
cal models and 

... to explore the semiclassical structure of the respec¬ 
tive quantized gauge field theory (also with dynam¬ 
ical fermions) 

within one and the same framework. 

In order to generate a large number of typical classical 
solutions of the Euclidean SU (3) Yang-Mills equations of 
motion for a given temperature and volume, Monte Carlo 
lattice gauge fields have been cooled with respect to the 
standard Wilson gauge action. The smoothed configu¬ 
rations have then be characterized with the help of an 
0{a^) improved field strength tensor 


F^^{x) 


lE 




( 10 ) 

where is the traceless and antihermitean part of 

the respective k x I rectangular Wilson loop. The sum is 
extended over the four loops of quadratic shape and each 
given size (n, n) in the /iJ^-plane which have the site 
X as a common corner. From the field strength one can 
evaluate the corresponding improved operators of action 
density s(x) and topological charge density qt{x) 


s[x) = ^tr(F^2^(x)) and (11) 

qt{x) = — e/ijypo-tr (L)ji/(x)E'po-(x)). (12) 


In our case of cooled gauge fields, the total topological 
charge Qt — deviates from integer values typi¬ 

cally by less than one percent. 

The starting configurations for our cooling search have 
been produced by Monte Carlo simulation at jd = &/= 


5.65. For TVt = 4 and 6 the Monte Carlo samples are 
in the confinement phase with T/Tdec ~ 0.91 and 0.6, 
respectively. In these cases the volume average of the 
Polyakov loop vanishes. Under the influence of cooling 
this average starts to diffuse. A large fraction of con¬ 
figurations will approximately remain in the vicinity of 
maximally non-trivial holonomy until each turns into a 
classical solution. Various lattice sizes 4 x 12^, 6 x 12^, 
4 X 20^, 6 X 20^, and 12^ have been chosen in order 
to study how the yield of classical configurations (pop¬ 
ulating different parts of the caloron parameter space) 
depends on the lattice size, in particular on the aspect 
ratio. Starting from the deconfined phase (T > Tdec), 
mainly fields with trivial holonomy have been obtained 
by cooling. Throughout this paper we restrict ourselves 
to the confinement phase since we are mainly interested 
in KvB solutions with generically non-trivial holonomy. 

To what extent a lattice gauge field can be regarded to 
be classical is measured by the violation of (anti) selfdu¬ 
ality 5p, defined by 


5f 


E 


^s(x) - |gt(a;)| 


(13) 


For a strictly selfdual or anti-selfdual gauge field this 
quantity vanishes. 

If the caloron constituents are well-separated the ac¬ 
tion density of the caloron becomes static in xq. An 
auxiliary observable which we call non-staticity 


W Ex |s(a^ + 6) - s{x) 
4 Ex'S(a^) 


(14) 


has been introduced to distinguish between dissociated 
monopole constituents (for which the action density is 
static) on one hand and the non-dissociated case (which 
gives rise to normal time-dependent instanton lumps) on 
the other. The factor iVt/4 compensates for finite dis¬ 
cretization as was chosen in previous work for SU (2) 
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For classical solutions, evaluating the continuum action 
density Eq. 0 on a grid of points, the following bench¬ 
marks have been calculated for a symmetric set-up of the 
constituents. Assuming that the holonomy is maximally 
non-trivial, i.e. tr('Poo) = 0, at least two constituents 
are discernible as maxima of the action density if 5t 
is lower than some bifurcation value 51 = 0.27. If the 
non-staticity is even smaller than 5^* = 0.06 , then it 
is possible to resolve all three constituents in the ana¬ 
lytically known action density profile of the SU{3) KvB 
caloron |4n|. 

Unlimited cooling with the standard Wilson action will 
always drive the field towards a topological trivial vac¬ 
uum IQtl ^ 0. In order not to loose all approximately 
classical fields the cooling process is monitored and a 
stopping condition is chosen. Using the non-staticity 5t 
and the violation of (anti) selfduality 6f we have em¬ 
ployed two different stopping criteria which allowed us 
to catch an abundance of (multi-)caloron solutions in or¬ 
der to explore the moduli space of calorons. Cooling is 
stopped, 

(A) if the violation of (anti) selfduality Sp passes 
through a minimum with the constraint Sp < 0.1,^ 

(B) if either 5p or the non-staticity 5t passes through a 
minimum, depending on what happens first. Addi¬ 
tionally the (less restrictive) constraint Sp < 0.2 is 
applied. 


Both criteria result in nearly (anti) selfdual gauge fields 
representing various topological sectors. In particular, in 
most cases configurations with some integer topological 
charge Qt and action S ^ S'inst \Qt\ (implying a slight 
violation 6p) are found. 

With the traditional choice (A) gauge fields are cooled 
until the equations of motion are fulfilled as good as pos¬ 
sible and only violated to a small extent. There are some 
reasons to argue that this stopping condition is too strict. 
Superpositions of semiclassical objects (e.g. instanton 
and anti-instanton) are no strict solutions of the equa¬ 
tions of motion. Nevertheless, a mixture of both must 
exist in a realistic model of the vacuum structure. On the 
other hand it is impossible to model the moduli space for 
an infinite volume. The numerical investigation can 
only be done for a discrete torus, where it is known that 
a \Qt\ = I field can exist only at the expense of violat¬ 
ing the (anti) selfduality |^. Accepting a certain viola¬ 
tion of the equation of motion constitutes a compromise 
between the moduli spaces of the torus and the finite 
temperature space-time x Furthermore lattice ar¬ 
tifacts of the action will bias the semiclassical content of 
fields. Hence one should not insist on perfect solutions 
of the (lattice) equations of motion. Stopping condition 


® This criterion is found to be practicall y equiv alent to a stopping 
condition used in our previous papers mui. 


(B) allows to stop in an earlier stage. In particular this 


results in an enrichment of static, approximate classical 
configurations. 


For a cooled gauge field configuration we determine 
the holonomy Voo as an average over a spatial region Va , 
where the action density is minimal, 


'Poo — 






diagonal 


xGVc. 


= diag(e*'^i(^\e*‘^^(“),e*'^^(=^)) (15) 

y r 


xGVc. 


with 

W 

Pix) = Yi U(^xo,s),o- ( 16 ) 

Xo = l 


The domain is taken as the fraction of a = 10 % of the 
3-dimensional volume for which the spatial action density 
s{x) (the XQ-average of s{x) with fixed x) is smaller than 
everywhere in the remaining region. To get a meaningful 
result for the Poo operator it is necessary to choose a 
suitable gauge. We decided to take the Polyakov loop 
diagonalizing gauge, i.e. we average over the eigenvalues 
of P. From the observed Polyakov loop operator P{x) 
not only the holonomy Poo (and the Hm’s) is determined, 
but also the (monopole) positions yn and consequently 
the overall number of constituent monopoles. 

Indeed, with the diagonalized Polyakov loop operator 
['P(a;)]j;„,gQjjg^j we identify the constituents by determin¬ 
ing the positions Y where two eigenvalues approach each 
other. For the caloron solution in the continuum they 
would exactly coincide. On the lattice this condition has 
to be slightly relaxed. Thus, in our numerical work we 
fixed monopole positions Y through the local minima of 
the function 


/(f) = min (17) 


in the 3-dimensional volume with the additional con¬ 
straint that some (j) exists such that 


min 

4>G[0,27t) 


L{Y) - (2e*‘^ -be“^*‘^)/3 


< 0.05. 


This ensures that (approaching some (f>) the eigenvalues 
of P are sufficiently degenerate at Y. To avoid spuri¬ 
ous monopole positions, only monopole constituents with 
s(U) > 0.1maX;js(f) are taken into account. 

For individual configurations we also compute the 
spectrum of the clover-improved Wilson Dirac operator, 
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example 1 - global observables 

S/Slnst 

1.14 

Qt 

-1.00 


-0.36 



0.00 



0.36 , 


n+ — n_ 

1-0 


CO 

1 

o 

CO 

t-H 

6 f 

0.14 


TABLE I: Properties of a Qt = — 1 caloron gauge field shown 
in Figs. El El and El 


defined as 

4 

y) — ^x,y ^ ^ ^ j^(f^4 0^/j,) Ux,(i H" 

(14 + 7/i) 

2 , 

lp{x,y) = 4{x,y) + -KCsv,cr^,,Fi_,^{x)Sa;,y 
with ^ [7 m>7i^] ■ 

The parameters are chosen to be k = 1/8 (massless) 
and Csw = 1 (tree-level improved). For all spatial direc¬ 
tions the boundary conditions of spinor fields are peri¬ 
odic, for the temporal direction one has variable ones: 
'I'z(2;o + b,x) = '^z{xo,x). Eigenmodes are com¬ 

puted using the ARPACK package j^. 


IV. PROPERTIES OF LATTICE CALORONS 


Under the conditions described above, different classi¬ 
cal fields with | Qt | = 1, • ■ ■, 7 embedded in an arbitrary 
holonomy are obtained from cooling. Since the Monte 
Carlo simulation is performed in the confined phase, the 
yield of (almost) trivial holonomy configurations is low. 
Here we present and explain examples for \Qt\ = 0,1,2. 
In order to illustrate these solutions, we provide tables 
with gluonic and fermionic observables. 

Our first example of a SU{3) classical lattice gauge 
field has topological charge Qt = —1 and non-trivial 
holonomy (cf. Tab.Pl. It has been generated by cooling 
with respect to the standard Wilson acti on o n a lattice of 
size 4 X 20^ with the stopping criterion (B) Due to the 
non-trivial holonomy the gauge field contains three (al¬ 
most equal) massive constituent monopoles, which can be 
identified by the monopole criteria (a-d). These criteria 
coincide since the monopoles are sufficiently separated. 
In Fig. 13 this is made visible in the form of isosurfaces 
of the action density, of the function f{x) according to 
Eq. 11711 showing points, where two eigenvalues of 'P{x) 
come close to each other, as well as isosurfaces of the 



FIG. 3: Caloron with non-trivial holonomy (example 1), topo¬ 
logical charge Qt = —1 and separated constituents. We show 
isosurface plots for the action density (u.L), for the function 
f{x) (Eq. inj) localizing the monopoles, for the fermionic 
density with boundary conditions x = 1/6,1/2, 5/6 

drawn with different colors/brightness at fixed xq (1.1.) and 
fixed Xi (l.r.). 


fermionic density for varying fermionic boundary con¬ 
ditions. Evidently, the significant features of the KvB 
caloron of this first example are visible in Fig. 13 but the 
field is not perfectly anti-selfdual (compare with the Sp 
valu e in T ab.P). This effect is due to the stopping condi¬ 
tion (B) which terminated cooling at the minimal value 
of 5t with the constraint 6f < 0.2. 

Further cooling, until stopping condition (A) applies, 
would squeeze the constituents together such that no sep¬ 
arate constituents would be visible in the action or topo¬ 
logical charge density. Still, the Polyakov loop V and also 
Abelian monopoles would indicate a substructure of the 
single massive lump. Fermionic observables would give 
only a slight hint to the substructure, since the fermionic 
zero-mode density only slightly moves around inside the 
lump and is “breathing” under changes of the temporal 
boundary condition of the Dirac operator. It “delocal¬ 
izes” to an algebraic decay at z = y,i. This situation is 
not visualized. 


The clover improved Dirac operator with the back¬ 
ground caloron field contains one zero-mode with posi¬ 
tive chirality. The localization, measured by the inverse 
participation ratio IPR = V , strongly 

depends on the boundary condition and signals a delocal¬ 
ization IPR ^ I where the jumping of the zero-mode is 
expected to occur, namely at z = 0] 1/3; 2/3 (cf. Fig. 0 
lower right). The maximal localization is expected at 
x= 1/6; 1/2; 5/6. 

Since the limit |£| —> oo is not accessible on a fi¬ 
nite lattice one might question the possibility to define 
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FIG. 4: Caloron (example 1): scatterplot of V{x) (for the 
points with high enough action density) in the complex plane 
(u.L); eigenvalues of 'P{x) versus spatial action density (eigen¬ 
values are sorted —0.5 < /n < A*2 < Ms 5; 0-5, as indicated by 
different points (colors/styles)) (u.r.); Im(A) (eigenvalue of 0) 
(1.1.) and inverse participation ratio of the lowest eigenmodes 
(hr.) versus the temporal fermionic boundary condition z. 


the holonomy on the lattice. Nevertheless, as the upper 
right part of Fig. 0] shows, the eigenvalues of V have 
very small variance over the ZD points where the spatial 
action s is low. This property of the solution makes our 
definition of 'Poo a posteriori plausible. Furthermore, 
the eigenvalues of V approach each other at points with 
larger values of the spatial action. Since one does not ob¬ 
serve an exact matching, the monopole criterium relying 
on Eq. 03 is numerically better defined. Confirming the 
analytic expectation, the action density is static as indi¬ 
cated by the small value of St and resides in separated 
constituents. 

Fitting the lattice field by analytic expressions of the 
517(3) KvB caloron solution is not feasible. The numer¬ 
ical reconstruction of the gauge links C/x,/x in terms of 
the vector potential Afj,{x) is too costly. Moreover, the 
“latticized” analytic solutions are not yet adapted to the 
periodic boundary conditions. But for gauge invariant 
quantities like s{x) and fE'l'kz(a:) the available ana¬ 
lytical expressions have been easily used for a fit. The 
zero-mode density is (exponentially) localized only at one 
respective constituent for almost all angles z, minimizing 
in this way effects of the spatial boundary condition. This 
is shown as an example in Fig. Eland turns out to provide 
a correct fit of the profile over two orders of magnitude. 

Finally let us have a look at the charge-one caloron 
in a maximally Abelian projection which is carried 
out similar to Ref. ^3- We maximize Fu{D,g) = 
Re ti {g{x)Ux,ij.9H^ + A) with respect to g € 



FIG. 5: Fit of the analytic expression of z{x) to the 
actual density for z = 'bl% \ full lines means measured; dashed 
lines correspond the fit. 



FIG. 6: Abelian monopoles within a Qt = —1 caloron with 
non-trivial holonomy and separated constituents (example 1); 
colored arrows show electric (= —magnetic) field components 
of Abelian projected C/(l) x U{1) field; the isosurface shows 
the action density of the Abelian projected field. 


SU{3) and D , where D^^ii is a C/(l)^ C SU{Z) matrix 
of the form ZI = diag(e*''^b e“h7i+72)) ^ Consider¬ 
ing the metric |17 —i7p=tr {U — D)\U — D) the gauge 
transformed link U is required to be as close as possible 
to the diagonal matrix D. Technically this is performed 
using the Cabibbo-Marinari method. The static elec¬ 
tric (equal to imagnetic part for an (anti)selfdual field) 
components of the U{\) x 17(1) field strength are com¬ 
puted from the improved field strength tensor Eq. cnj, 
where the SU{3) links Ux,ti are replaced by the pro¬ 
jected Dx /i’s. The resulting field pattern is shown in 
Fig. El The maximally Abelian monopoles (sources of 
the electric and magnetic fields) reside inside the spheres 
where the Abelian field becomes singular. 

The second example of a selfdual lattice gauge field 
has Qt = 2 and also possesses a non-trivial holonomy 
(cf. Tab.ini). It is obtained by cooling with respect to 
Wilson gauge action and by applying stopping condition 


(B) Again the expected full number of 3 |Qt| = 6 



















example 2 - global observables 

S/Slnst 

2.04 

Qt 

2.00 


-0.35 


hi j 

-0.01 



0.36 , 


n+ — n_ 

0-2 


CO 

1 

0 

6f 

0.04 


TABLE II: Properties of a Qt = 2 caloron gauge field shown 
in Figs. 0and|^ 





FIG. 7: Caloron (example 2) with non-trivial holonomy 
and IQtl = 2; six constituents are visible in all plots: ac¬ 
tion density (u.L), monopoles by degenerate eigenvalues of 
P{x) (u.r.); for 2 : = 1/6, 1/2,5/6 (different col¬ 

ors/brightness) (1.1. and hr.). 



monopoles can be identified by the monopole criteria (b- 
d) - note that there are only 5 separate massive lumps 
visible as separate maxima in the isosurface plot of the 
action density and an application of criterion (a) is more 
subjective The decomposition into six monopole con¬ 
stituents can be seen in Fig. |7| 

By analyzing example 2 the results for |Qt| = 1 are 
confirmed and generalized. Hence lattice calorons, espe¬ 
cially those with non-trivial holonomy, show the typical 
properties of KvB calorons. Further cooling reduces the 
violation of the equation of motion and does not influence 
the constituent positions too much. However, there is a 


^ A more sophisticated cluster algorithm for example, would be 
able to find six clusters of topological charge 




FIG. 8: Caloron (example 2); scatterplot of Vix) for the 
points with high enough action density) in the complex plane 
(u.L); eigenvalues of V{x) versus spatial action density (u.r.); 
Im(A) (eigenvalue of 0) (1.1.) and inverse participation ratio 
of the two zero-modes (hr.) versus the temporal fermionic 
boundary condition 2 . 


striking difference to the case \Qt \ = 1 which is a con¬ 
sequence of the Nahm duality [^, which points out that 
classical fields with \Qt\ = 1 are unstable on the torus. 
In agreement with the Atiyah-Singer index theorem the 
Qt = 2 field supports two zero-modes with negative 
chirality. The 2 -dependence of the inverse participation 
ratio IPR of (x) and the low lying eigenvalues of 

the Dirac operator are shown in the lower part of Fig. 
The behavior of both zero-modes - the choice of basis is 
somewhat ambiguous - is similar to the |Qt| = 1 case. 
The Polyakov loop through the constituents aligns as it 
is predicted by Eq. O- The slight deviation, which can 
be observed in the upper left part of Fig. |51 might be 
due to finite volume effects. However, the upper right 
plot of the same figure shows that the holonomy Voo is 
again well defined (witnessed by the small variance of the 
eigenphases over the “asymptotic” subvolume V/t)- 
The third and last example of an almost classical gauge 
field has Qt = 0 but needs to be interpreted in terms of a 
caloron-anticaloron superposition. The holonomy is not 
maximally non-trivial, giving rise to a light, a moderately 
large and a heavy monopole constituent (mass fractions 
P « (0.1; 0.35; 0.55) cf. Tab. Ell). This 4 x 123 
field has been cooled with respect to the standard Wilson 
actio n and again obtained under the stopping condition 
(B) Only 4 < 3 |5'/5'inst| = 6 constituents can be 


identified as massive lumps of action (a), which is not 
surprising keeping the monopole masses in mind. 

The upper panels of Fig. E show the isosurfaces of action 
density s{x) and the topological charge density q{x). 
The interpretation in terms of a caloron-anticaloron su- 








































9 


example 3 - global observables 

S/S'in.t 

2.00 

Qt 

0.00 


-0.18 


hi { 

-0.08 



0.26 , 


n+ — n_ 

0-0 


6.0 • 10-2 

5f 

0.18 


TABLE III: Properties of a caloron-anticaloron gauge field 
shown in Figs. l^andlTiil 



FIG. 9: Caloron-anticaloron superposition (example 3) with 
S = 2S'inst , topological charge Qt = 0 ; action density (u.L), 
topological charge density (red=positive, blue=negative) 
(u.r.); TJ 75 \I'z(x) for z = 0,0.5 (in different col¬ 

ors/brightness for each would-be zero-mode), shown in 3D 
space (1.1.) and depending on Euclidean time xq (l.r.). 


perposition is supported by the relatively large value of 
Sp = 0.18 compared with example 2. Our choice of 
Sp < 0.2 barely allowed us to capture this approxi¬ 
mate solution when the non-staticity St was passing 
a minimum. Such fields can become ultimately stable 
( i5f —> 0 ) only if the selfdual and anti-selfdual centers 
are infinitely separated. Any finite separation results in 
an attractive force between the caloron and anticaloron 
that renders the superposition unstable. 

There is no zero-mode present in this case. Instead 
we consider the pair of non-zero modes (with the small¬ 
est non-vanishing imaginary part of the eigenvalue A). 
For this pair of modes one needs to have = 

0 . Hence, both eigenmodes, with ImA positive and 
negative, need to be localized simultaneously on the 
caloron and the anticaloron for a chosen boundary con- 




FIG. 10: Caloron-anticaloron superposition (example 3): 
scatterplot of Vix) (for the points with high enough action 
density) in the complex plane (u.L); eigenvalues of Vix) ver¬ 
sus spatial action density (u.r.); Im(A) (eigenvalue of 0 ) 
(1.1.) and inverse participation ratio of the lowest eigenmodes 
(l.r.) versus the temporal fermionic boundary condition z. 


dition. The characteristic jumping, presented for |Qt| = 
1,2 , can be seen for the positive or negative part of 
(x) . Both parts enclose lumps of action with 
the same sign of topological density. The flow of the in¬ 
verse participation ratio IPR of the eigenmode density 
in Fig. cni still shows two similarly localized eigenmodes 
(cf. circles and diamonds in the lower right panel). In 
the same plot it is visible that the modes only delocalize 
twice, which would support an interpretation in terms of 
an embedded SU(2) caloron-anticaloron. This is actu¬ 
ally not the case, e.g. the scatterplot of the Polyakov loop 
clearly exhibits a scattering in the full complex plane, 
which is not possible for a trivial SU (2) embedding. The 
interpretation in terms of a KvB caloron is still valid, 
however, the behavior of a superposition of fields with op¬ 
posite duality® differs from a fully selfdual or anti-selfdual 
caloron. 


V. ENSEMBLES OF LATTICE CALORONS 

To underline the importance of classical solutions with 
non-trivial holonomy we show on the left in Fig. llll a dis¬ 
tribution of the modulus of the average of L{x) , which 
would be equal to one for trivial holonomy Voc G Z(3) 
and smaller for non-trivial holonomy. It shows primar¬ 
ily that we are investigating properties of classical fields 


® Only locally the field is mainly either selfdual or antiselfdual. 
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topological charge 


FIG. 11: Histogram of |trPoo|/3 (left) and Qt (right) for 4 x 12® (stopping condition |(B) 


with a non-trivial holonomy and that these fields can be 
produced by cooling on the lattice. The distribution of 
topological charges on the right in Fig.^Jshows the avail¬ 
able range of topological sectors. The case Qt = 0 is 
excluded, because configurations of this class mostly pass 
cooling without being stopped before they reach very low 
action (the third example above is a fortunate exception 
just tolerated by 5f < 0.2). Both distributions only 
weakly depend on the stopping condition. The distri¬ 
butions of the non-staticity 5t and of the violation of 
selfduality {i.e. of the equation of motion) 5f naturally 
dep end o n the stopping condition. According to condi¬ 
tion (A) cooling proceeds until the equation of motion 
is maximally fulfilled. In this (actually late) stage the 
constituents are already squeezed. In the correspond¬ 
ing ensemble of cooled configurations one cannot find 
static calorons with separated constituents, at least for 
Qt| = 1 where squeezing is found to be a particularly 
strong effect. With the help of the non-staticity this role 
of the stopping criteria is visualized in the upper panels of 
Fig- El “ with stopping condition |(A)| (upper right plot) 
no static calorons with |Qt| = 1 < g; 1 ) are found, 

whereas the combined stopping condition (B) (upper left 
plot) gives both static (peak (a)) and non-static calorons 
(peak (b)) with |Qt| = I. This effect is much weaker for 
IQtl =2,...,4. 


For the general understanding it is also interesting to 
know how the violation of selfduality is distributed. 
For both stopping criteria this is shown in the lower pan¬ 
els of Fig. ^ Clearly, the case |Qt| = 1 appears ex¬ 
ceptional also in this respect since it has a positive lower 
bound for the violation of (anti)selfduality of « 0.02 in 
contrast to classical configurations with |Qt| > I. 


The previous figures describe the composition of the 
ensembles. In order to characterize better the composi¬ 
tion of the configurations one should choose observables 
closer related to the KvB caloron solution. It would be 
most natural to count the constituent monopoles inside 
the caloron configurations. For this purpose we apply the 
definition of a monopole as sitting at Y by two coinciding 




FIG. 12: Distribution of non-staticity St (upper) and of the 
viola tion of (ant i) sel fduality Sf (lower) for stopping condition 
|(A)| (right) and |(B)| (left). 


eigenvalues of V{Y), which was explained in Chapter |lIII 
Additionally we show in Fig. 11 .‘II the number of domi¬ 
nant {i.e. |9t(F")| > 0.1 maxa; |gt(a;)| ) local extrema in 

the topological charge d ensit y for the ensemble produced 
with stopping condition (B) 


Both methods, the number of monopoles provisionally 
defined as local maxima of the action/topological charge 
density or by our actually adopted monopole definition, 
expose (cf. in Fig. I13II always more than |Qt| lumps, a 
situation that would be expected for a caloron with mass¬ 
less or close-by constituents. With the actual monopole 
definition, on average 3|Qt| monopoles are found in each 
topological sector. This is independent of the cut, which 
can be applied to the ensemble to make the monopole 
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1 2 3 41 2 3 41 2 3 4 

S/S. . 


FIG. 13: Number of constituents for 4 x 12® ensemble with standard deviation (bars); full squares by “coinciding of eigenvalues” 
criterion, open circles from local maxima in the 4D action density; (left) full ensemble; (middle) with L^o < 0.2 and (right) 
with St <0.1 and Loo < 0.2. 


counting better defined®. By applying more strict cuts 
to the ensemble (with respect to the holonomy and non- 
staticity, from left to right in Fig. ESI) the standard de¬ 
viation of the number of monopoles for given S/S-mst 
decreases while on average the expectation to find 3 |Qt | 
monopoles is reconfirmed. 


The splitting into static and non-static configurations 
in the distribution of the non-staticity (cf. upper left 
in Fig. PHl for stopping condition |(B)| is also reproduced 
looking at the distribution of the average monopole dis¬ 
tance in Fig. d In this Fig. the sample is restricted 
to configurations with \Qt \ = 1 where 3 monopoles have 
been observed. The correspondence between the clus¬ 
tering into static and non-static sets of configurations 
in Fig. IT^ and Fig. is highlighted by labelling those 
clusters of events (in both figures) as (a) and (b). They 
have emerged from the cooling process either because of 
the non-staticity went through a minimum (cluster (a)) 
or the violation of the equation of motion was minimal 
(cluster (b)). 

In Fig. dlit becomes clear by comparison that the first 
cluster (a) is well represented by calorons with con¬ 
stituents uniformly distributed on the lattice. The aver¬ 
age distance is comparably large and consequently they 
are static. On the other hand, configurations forming 
cluster (b) have emerged from cooling because of a min¬ 
imal violation of the equation of motion achieved in a 
later stage of cooling. The constituents are squeezed and 
the action density has become non-static, which can be 
monitored during cooling by increasing values of St- 
The correlation between non-staticity and monopole 
distance qualitatively satisfies the analytical expectation 
while the effect of the stopping condition reflects the ex¬ 
pected, unavoidable side-effect of cooling with respect to 


the standard Wilson action. 

To understand what changes when the temperature is 
lowered at fixed volume and approaches zero tempera¬ 
ture (on the symmetric torus) we have studied the same 
properties for an ensemble on a 6 x 12® and a 12^ lattice. 
Nr and together (by the aspect ratio of the lattice) 
characterize to what extent finite T in infinite volume, 
S'® X R®, is modeled. Only if TVo- 3> Nr the caloron finite 
volume effects (caused by overlapping constituents) can 
be neglected. By sending Nr —> iV^ and keeping the 
size (and distance) of the object small compared to the 
box length, more instanton-like classical solution are 
expected. 


In our earlier SU (2) lattice studies a recombination 
of dissociated caloron constituents into a single non- 
dissociated object, still with non-trivial holonomy, has 
been observed in the limit Nr N^ [3 We repeat 
this study for SU{3) by using lattices with aspect ratios 
Nr/Na = 6/12 and Nr/N^ = 12/12 in creating ensem- 
bles of configurations cooled until the stopping condition 


(B) is satisfied. In complete analogy to the SU{2) case 


static calorons disappear if the aspect ratio Nr/N^- is 
increased (upper part of Fig. ITKll . In Fig. ^Ithe num¬ 
ber of constituents, found for different topological sectors 
IQtl , for Nr = 6 is shown. Different cuts with respect 
to the holonomy \Voo\ and non-staticity St are applied. 
Despite the bigger aspect ratio 3 IQtl monopoles are still 
visible whereas a smaller number of local extrema is ob¬ 
served in the topological charge density, closer to one 
lump per unit charge, if no further cuts with respect to 
holonomy and non-staticity are applied. Since the inter¬ 
pretation of the Polyakov loop (corresponding to a choice 
of a “time” direction) would be ambiguous, this has not 


Without this cut one finds a large number of spurious monopoles, 
in particular if the holonomy is not sufficiently non-trivial. 


^ A similar observation has been made employing the me thod of 
“adiabatic cooling” which uses an anisotropic lattice l38l . 
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average monopole distance 


FIG. 14: Distribution of the average pairwise monopole distances compared with the result of an uniform distribution of 3D 
points (left). Average monopole distance versus the non-staticity <5*. 



FIG. 15: Distribution of non-staticity (upper) and topological 
charges (lower) for the 6x 12^ (left) and 12"^ (right) ensembles. 


been checked on the symmetric lattice. We would like 
to remind the reader that Gattringer and Pullirsch |4^ 
have found on the symmetric torus, similar to the ob¬ 
servations at T 0 [S^, jumping of the zero-mode for 
Qt = ±1 Monte-Carlo configurations under a change of 
boundary conditions. In the light of the present studies 
we would doubt an interpretation in terms of jumping 
between fractionally-charged constituents. 


VI. CONCLUSIONS 

We have applied cooling to SU{3) Monte Carlo helds 
in the confined phase at finite temperature. The typ¬ 
ical KvB behavior of these helds was shown for exam¬ 
ples with \Qt\ = 1,2 and Qt = 0 by analyzing held 


strength, Polyakov loop, action and topological charge 
prohles and the behavior of low-lying modes of the Dirac 
operator. This part was complemented by studying the 
content of the whole ensemble, especially with respect 
to the monopole content of solutions in each topological 
charge sector. The outcome is biased by the choice of 
lattice action - we have used the standard Wilson action 
- which can be compensated by a suitable stopping con¬ 
dition for the cooling iteration. Still, \Qt \ = 1 plays an 
exceptional role. However, in agreement with the analyt¬ 
ical expectation we hnd 3 jQtl constituent monopoles 
(for |Qt| = ) which are dehned by the coinci¬ 

dence of two eigenvalues of V{x) . 

Zero modes of the Dirac operator with adjustable 
boundary conditions are an useful filter to explore the 
semiclassical structure. Chirally improved Dirac opera¬ 
tors are necessary for the study of its low-lying modes to 
be applicable directly to unsmeared Monte Carlo conhg- 
urations. Using the full feature of exploring the spectral 
flow depending on the boundary condition leads to high 
computational costs such that only small samples of con¬ 
figurations can be studied. 

Uncovering the topological structure requires one to 
introduce notion of smoothness. Therefore smearing has 
been applied to explore the topological density or the 
action density. Both observables can give valuable infor¬ 
mation about the lattice configuration already after some 
amount of cooling or smearing. Nevertheless, one might 
miss significant features of KvB calorons by relying only 
on those. Having the connection between the holonomy 
and deconfinement-confinement phase transition in mind, 
the Polyakov loop operator is another very important and 
valuable tool to investigate the classical properties. The 
signature of a monopole is that L{x) = (1/3) tiV{x) is 
close to the envelope of the range of the Polyakov loop 
(see the upper left panel of Fig. El). But it is unlikely that 
this signature can be directly observed in an early stage 
of smearing or cooling. 

In future investigations we will use these tools to ex¬ 
plore the properties of finite temperature lattice field 
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no cut imposed 



|L I <0.2 |LJ<0.2 and 5^ <0.2 



FIG. 16: Dissociation into constitnents from local maxima of action and monopole counting (defined by coinciding eigenvalues) 
for 6 X 12® lattice with different cuts imposed. 


theory near the SU (3) deconfinement-confinement phase 
transition by applying different smearing techniques. 
Lattice studies can hopefully help to identify the disso¬ 
ciation of calorons with non-trivial holonomy if the tem¬ 
perature is lowered below the critical one, in other words, 
to verify the semiclassical results. 
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APPENDIX A: AUXILIARY FIELDS 

The auxiliary fields for the construction of SU{3) 
|Qt| = 1 calorons are [^: 

i>{xQ,x) = ]^iv[AzA2Ai)-cos{2TTXo), 

Pm+l I 

'^m+1 




using the notation chm = cosh(27rr'm'rm), shm = 
Smh{2T:Vrnrm), Tm = \x-ym\ and Pm = Vm - Vm-l- 
The fermionic density is derivable from the Green’s 
function 


fx{z,z') = 


^27Tit{z — Z^) - 


e sinh{27Trm{z — z)) 


rmi^ 

+ {Vm{z')\TT{A3A2Al)\Wm{z)) 


with 


= -wl^iz) = sinh(27rrm(z - pm)), 

Vmiz) = w^(z) = cosh(27rrm(z - ^m)). 

The construction of the caloron field requires the 
knowledge of the following matrix helds Ca and (j)- 

iC(x)jrik — CmVap Ck Pk\ 

) mk ~ Smk ~ CmCk fx{Prm Pk)- 


The anti-selfdual ’t Hooft symbols are fjap = d[Q,cr^]- 
The permutation tc{A 3 A 2 Ai) is A 3 A 2 A 1 , A 1 A 3 A 2 
or A 2 A 1 A 3 for m = 1,2 or 3, respectively. In the 
expression for / the assumption pm <z'<z< Pm+i 
has been made. This can be generalized by fx{z',z) = 
f*{z,z'). If the constituent positions ym lie in the Xi — 
X 2 plane, one can use Cm = (jp^MPm “ Pln)l \/27r|pm|. 
In general one has pm = —T^CmfCm with the Pauli vector 
T = (ri,T 2 ,T 3 ). The temporal shift yo is set to zero 
and b = 1/T = 1 throughout these formulas. The 
angles pm are determined by the holonomy Voo — 
exp(27ri diag{pi, p 2 , P 3 )). 






